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Abstract 



We report on a Hartree-Fock approximation study of the meron 
pseudospin-texture excitations of the broken-symmetry incompressible ground 
states of double-layer quantum Hall states at v = 1. We have obtained results 
for meron core energies and core sizes. We use these to estimate the charge 
gap which determines the activation energy for dissipation in the quantum 
Hall effect. 
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Recent work from several different points of view has led to the realization that 
double-layer two-dimensional electron systems (2DES) in a strong magnetic field can, under 
appropriate circumstances, exhibit an unusual broken symmetry in which interlayer phase 
coherence develops in the absence of interlayer tunneling. This broken symmetry has useful 
analogies with the spin-ferromagnetism which occursi in a single-layer 2DES in the quantum 
Hall regime, for example, at odd integral Landau level filling factors, v. In the pseudospin 
representation commonly used for the (spin-polarized) double-layer system, spontaneous in- 
terlayer phase coherence is equivalent to the spontaneously broken U(l) symmetry of an 
easy-plane ferromagnet.ii In the limit of zero separation d between the layers, the (spin- 
polarized) double-layer system with no tunneling is equivalent to the single-layer system with 
no Zeeman coupling; in this limit the ferromagnetism of the pseudospin representation for the 
double-layer 2DES becomes isotropic. Interesting physical consequences follow from sponta- 
neous interlayer phase coherence, including the occurrence of a finite temperature Kosterlitz- 
Thouless phase transitionBBi, superfluid-like behavior in the pseudospin channel,@@i and a 
commensurate-incommensurate phase transition driven by an in-plane component of the ex- 
ternal magnetic field in the presence of weak interlayer tunneling.0B0 Argueably, the most 
remarkable property associated with the broken symmetry is the occurrence of highly collec- 
tive, topologically nontrivial pseudospin texture excitations which carry a physical charge 
equal (in units of the electron charge) to the topological charge of the texture. 0i These 
charged pseudospin textures are analogous to the skyrmion charged spin textures of single- 
layer systems BOB whose presence has recently been detected experimentally.lll In this paper 
we report on a Hartree-Fock calculation which estimates the energy of these excitations for 
double-layer 2DES's at v = 1. 

The fractional quantum Hall effect, an anomaly in the transport properties of 2DES's 
in a strong magnetic occurs when the system has a charge gap, i.e., a discontinuity 

in the chemical potential, at a critical density which depends on magnetic field. (It is the 
magnetic field dependence of the critical density which leads to the charge carried by the 
spin and pseudospin textures.!) At the critical density the dissipative conductivity, a xx , 



has an activated temperature dependence with an activation energy which, for the case of 
interest here, is half the charge gap. In single-layer systems, it is knownllll! that the lowest 
energy charged excitations at v — 1 (but not at larger odd integral values of v) are skyrmion 
spin-textures. It follows from continuity that in spin-polarized double-layer systems the 
charge gap must also be given by the energy of its charged spin textures, at least if the layer 
separation, d, is not too large. In the double-layer case, where the broken symmetry is planar, 
vortices might be expected to play an important role. Indeed, it has been proposed^! that 
the charged pseudospin textures of double-layer systems consist of meron-antimeron pairs. 
Merons are vortex-like pseudospin textures where the pseudospin tilts out of the easy plane 
in the vortex core; neither the core size nor the core energy are fixed by the parameters which 
characterize the system at long length scales!!! and both must be determined by microscopic 
calculations. Merons have an infinite energy and carryHi charge ±e/2, but meron pairs with 
like charges and opposite vorticities carry charge ±e and have a finite energy. In this paper 
we evaluate the meron-core energy and size using a Hartree-Fock approximation and use the 
results to estimate the energy of the charge ±e excitations of double-layer systems and to 
limit the range of validity of the meron pair description of these excitations. 

In a recent paper,i Moon et al. introduced single-Slater-determinant microscopic wave 
functions for both skyrmion and meron states near v = 1, for which the spin and pseudospin 
textures have the appropriate topological charges. The form of these wave functions provides 
an understanding of the connection between topological charge and physical charge from 
a microscopic physics point of view, as we discuss below, but the wavefunctions were not 
energetically optimized. For the single-layer case, Fertig et aZ.EH were able to find the optimal 
single-Slater-determinant wavefunctions for a skyrmion spin texture as a function of the 
strength of the Zeeman coupling by solving a set of Hartree-Fock equations. In this paper 
we generalize the Hartree-Fock calculations to the case of double-layer systems and find the 
optimal single-Slater-determinant approximation to the meron state as a function of d. The 
technical details of our calculation are very similar to those of Fertig et al. However, meron 
states are quite different from skyrmion states and have richer physics as we demonstrate 
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below. 

We assume that the Landau level spacing is so large that all electron stay in the lowest 
Landau level (LLL), hence we only need to keep one body states in the LLL. To take advan- 
tage of the circular symmetry of the meron state, we use the symmetric gauge in which the 
one body wave functions in the LLL have the form 4> m (r) = (2Ti2 m m\)~ 1 ^ 2 z m exp(— 12;| 2 /4), 
where z = x + iy is the complex coordinate, m is the angular momentum quantum number 
and lengths are expressed in units of the magnetic length. [I = hc/eB where B is the 
magnetic field strength). In the limit d = (which is equivalent to the real spin case), inter- 
and intralayer Coulomb interaction have the same form and the system is invariant under 
pseudospin rotation.!! In the Hartree-Fock approximation, the ground stateB'i at v = 1 in 
the standard pseudospin representation is completely polarized along an arbitrary direction 
in the x — y plane: 

l*o> = II(4f l T + ^i|)|0>- (1) 

Here |0) is the vacuum state, c^,^ creates an electron in the m state in the upper (lower) 
layer, N is the number of electrons, and <fi is an arbitrary phase. (In the following we choose 
= so that the pseudospin is oriented in the x direction in the ground state.) Note that 
the orbitals are peaked farther from the origin as the angular momentum increases. For 
d — 0, the state in Eq. ([!]) is the exact ground state! and we expect it to remain accurate 
as long as d is not too large. 

The single-Slater-determinant states for merons centered at the origin with vorticity +1 
and with charge — e/2 and +e/2 have the formi: 

N-l 

l*+l-|> = II ( MmC m| + Vm4+U)l°>> ( 2 ) 
m=0 

N-2 

= C ll II (-^ C lnT + U mC ] m+H )\®). (3) 
m=0 

Here u m and v m are variational parameters which, for iV — ► oo, have the same values in both 
wavefunctions because of the particle-hole symmetry of the HamiltonianEI, and which satisfy 



the constraint \u m \ 2 + \v m \ 2 = 1 so that the wavefunction is normalized. The easy-plane 
anisotropy of the double-layer system requires that 

lim W = 1, (4) 

| u m | 

so that at the charge density is the same for both layers far enough from the meron center 
and the electrostatic energy cost of the meron texture is finite. @B The charges of the two 
states are ±e/2 as can be verified by evaluating the differences in mean occupation numbers 
between the meron states and the ground state. 

In this paper we assume, without loss of generality, that u m and v m are real. Far from the 
origin, single-particle orbitals with adjacent angular momenta have nearly identical magni- 
tudes near their peaks where they satisfy m+1 (r)/0 m (r) ~ exp(i0) where is the angular 
coordinate. It follows that for large m, the single-particle state in the Slater determinant is 
effectively proportional to the pseudospinor whose azimuthal angle is equal to the angular 
coordinate: 



m = 7=2 



( cos(fl/2) ^ 



(5) 



i sin(0/2) exp(i(J)) j 

(We have imposed the normalization constraint by setting u m = cos(0 m /2) and v m = 
sin(# m /2). 9 is the polar angle for the pseudospin contribution from orbital m.) At large 
m, 9 = 7r/2, the pseudospinor is confined to the x — y plane and has vorticity +1. At the 
center of the meron only m= o(r) ^ so that the pseudospinor points up for the charge 
— e/2 meron and down for the charge e/2 meron. Meron states with opposite pseudospin 
vorticity can be generated by interchanging the roles of the two layers. 

The Hartree-Fock equations we solve determine the u m and v m which minimize the 
expectation value of the Hamiltonian in states of the form of Eq. (|2|) and Eq. (|3]). We 
perform our calculation at finite iV (up to N = 60) and to mitigate finite-size effects we 
include in the Hamiltonian the external potential from a fixed background which neutralizes 
the ground state charge density. The full Hamiltonian of the system is 
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H = - V V V aia2 c f c f c c 

rnim2mzrrn a\02 

— , N ~ 1 1 

~ ^ y ^mrr^ma \ ^ ] 7y {V mnmn ~\~ V mnmn )} , (6) 

m cr n=0 

where a, =f or j; = = V A is the intralayer interaction and = is the 
interlayer interaction. For the sake of definiteness we have taken the electron layers to have 
negligible thickness so that 

|ri - r 2 | 



/g2 
dridr 2 (j)* mi (ri)^ 2 (r 2 ) r0 ma (ri)0 m4 (r 2 ) , (7) 
ri — ro 



^Snwrnsfw = y dr i dr ad( r i)C 2 ( r 2)-r= 0m 3 (ri)0m 4 (r2). (8) 



and 



ki - r 2 | 2 + d 2 

The one body term in (|]) is due to the interaction of the electrons with the neutralizing 
background charge; it reduces to an irrelevant constant for N —>■ oo and, although we include 
it in our numerical calculations, we omit it to free the discussion below from unnecessary 
clutter. 

An elementary but somewhat tedious calculation gives the following result for the ex- 
pectation value of the Hamiltonian in the charge — e/2 meron state: 

Y N-l 



(^+1,— i I -^1 ^+1,— \ ) r\ ^ y \^vnV J n^y^mnmn Kraifim) 

m,ra=0 

+v 2 v 2 (V A -V A \ 

1 u m u n \ y m+l,n+l,m+l,n+l v m+l,n-\-l,n+l,m+l) 

+2u m v n V mtn+lrn ^ n+1 — 2u m v m u n v n V m ^ n+1 ^ n ^ m+ i\. (9) 

Minimizing Eq. (J|) with respect to 9 m gives the following equation which must be solved 
self-consistently for the pseudospinor polar angles, 

tan(0 m ) = A m /B m , (10) 

where 

2V-1 

A m = -2 Kf,n + l,n,m + l ^n) (H) 
n=0 



and 

JV-1 

B - \V A — V A — V A 4- V A 

1J m / y [ v mnmn mnnm m+l,n+l,m+l,n+l ' m+l,n+l,n+l,m+l 
n=0 

' m,n+l,m,n+l n,m+l,n,m+l V ^/v mnmn m+l,n+l,m+l,7t+l 

_T/ A - - I/ £ - 1/ B "ll (121 

v mnnm ' m+l,n+l,n+l,m+l m,n+l,m,n+l r n,m+l,n,m+l/J • \ / 

Similar equations can be obtained for merons with charge +e/2. The various terms in A m 
and _B m can be understood rather simply. For large m the matrix elements are insensitive to 
unit shifts in the angular momentum indices so that the constant terms in the denominator 
vanish rather quickly. These unit angular momentum shifts are important in determining 
the the structure of the meron core but we will temporarily ignore them for the following 
qualitative discussion. The term proportional to cos(# n ) in the denominator is the sum 
of intralayer exchange matrix elements and the difference of direct matrix elements for 
intralayer and interlayer interactions. If V A were equal to V E Eq. (|10|) would be satisfied for 
any constant 9 n , corresponding to the pseudospin isotropy at d — 0. The difference between 
intralayer and interlayer matrix elements which appears in the denominator of Eq. ( [TOD 
gradually forces 6 m to 7r/2 as we move away from the meron core. 

The most important result of our calculation is shown in Fig. where we plot our 
estimates of the chemical potentials at densities just larger and just smaller than the critical 
density at which the incompressible state occurs, fi^p and fiJ, IP respectively. These chemical 
potentials were obtained from the optimized meron state energies by the following procedure. 
For each value of d the difference between the meron state energy and the ground state 
energy, was evaluated at a series of iV values. The meron core energies, E^ c were 
extracted by fitting our results to the expected form 

E± = E± c + Trp E HR/Rmc), (13) 



where R = y 2(iV — 1)£ is the radius of the system, is the Hartree-Fock approximation 
for the in-plane pseudospin stiffness!, 

^ = ^2 l°° dkk*V E {k)e- h2 '\ (14) 
62tc z Jo 



Ve(&) is the Fourier transform of the interlayer interaction. The meron core radius 
R mc (which is identical for merons with charges ±e/2) was defined by m z (R mc ) = 
ATr£ 2 (S z (R mc )) = 0.1, where S z (r) is the z component of the pseudo spin density operator.! 
The energy of an individual meron is infinite because of the large distance contribution to 
the gradient energy, but the energy of a meron pair of opposite vorticity is finite. The 
optimal separation!, R*, for the meron pair is determined by minimizing the sum of the 
Coulomb energy and the gradient energy. For R >> R mc we findi that R* = e 2 /8ttpe and 
E MP = 2E mc + 2ttpe(1 + \n(R* / R mc )) . In Fig. we plot R* and R mc as a function of 
d; the meron pair picture of the charge excitations is valid only when R > Rmc which is 
satisfied for d/£ > 0.6. Results at smaller values of d/£ in Fig. |J and Fig. are plotted 
as dotted lines. E MP is the energy to make a meron pair at fixed total electron number by 
adding or removing charge from the edge of the system to compensate for the charge of a 
meron; to add or remove an electron it is necessary to change the electron number in the 
incompressible state. It follows&H that 

Vmp = ±E MP + e {d) (15) 

where eo(d) is the energy per electron in the incompressible state. We remark that as 
a consequence of a particle-hole symmetry which applies in the thermodynamic limit 
I 1 mp + I 1 mp = 2eo(0). We have used this identity to check the reliability of our numeri- 
cal procedures, particularly those involved in the extraction of the meron core energy from 
the size dependence of the meron energy. For dji >~ 1.2 we do not believe that our esti- 
mates are reliable; results in this range of dji are also shown as dotted lines in Fig. [1| and 
Fig. JJ. The problem at these values of dji is connected with the fact that the Hartree-Fock 
ground state at a pseudospin-density-waveEl for dji >^ 1.2. Within the range of 

d that the meron pair picture of charge excitations is valid, the meron-pair states are clearly 
energetically more favorable than naive Hartree-Fock single particle states (see Fig. and 
our results agree with exact diagonalization studies^ qualitatively. 
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are only interested in meron excitations in the broken U(l) symmetry phase, where the 
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FIGURES 

FIG. 1. Estimates of the chemical potentials (in units of e 2 /(e£)) just above (Mmp) an d below 
(fJ-Mp) ^ = 1 as a function of d/£, assuming that the meron pair state is a good description of 
the charged excitations. The dotted portions of the lines are for the ranges of d were the meron 
pair picture is invalid and our estimate is unreliable. The difference between fi^p an d A*mp * s 
the charge gap of the system at v = 1. For comparison we also plot the chemical potentials for 
the Hartree-Fock single particle excitations (dashed lines) which cost more energy than the meron 
pairs. 

FIG. 2. Estimates of the optimal meron core size R mc and optimal separation R* between 
merons in a charged meron pair state as a function of dji. The dotted portions of these curves are 
in ranges of layer separation where the meron pair picture is not valid. 
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